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This paper is the result of setting up GRENOUILLE in the Nonlinear Dynamics Laboratory at 
the University of Maryland at College Park. With the experience acquired in the process of setting 
("/") , up GRENOUILLE, this manual was compiled from literature and from hand-on experience to serve 

as a quick guide, a step-by-step help to construct GRENOUILLE and to understand some of its 
' basic principles. 
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Q ' I. INITIAL CONSIDERATIONS 

(N ' 

Frequency resolved optical gating; (FROG) consists of an autocorrelator with its output going into a spectrometer 

, | connected to a CCD camera pi S 0, 0] . The CCD camera records a spectrogram called the FROG trace (FigO- 3 

I dimensional graph with intensity as a function of time delay on the horizontal axis and wavelength on the vertical 
axis). From this spectrogram, the full electric held (intensity and phase) can be reconstructed with the aid of a 
numerical interactive algorithm. The mathematical form of the FROG trace depends on the nonlinearity used to 
q generate it. For a second harmonic generating (SHG) crystal, the FROG trace is proportional to Equation where 
^ ' the proportionality constant is divided out when the results are presented. 

O ■ 

Ifrog{u,t) = \ dtE(t)E(t-T)exp(iut)\ 2 (1) 
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GRENOUILLE (grading-eliminated no-nonsense observation of ultrafast incident laser light e-fields) is the simplest 
SHG FROG device ever built p|. A typical FROG device has a delay line, the nonlinear medium and a spectrometer. 
GRENOUILLE replaces the delay line by a Fresnel biprism and combines the spectrometer and the nonlinear medium 
. into a thick SHG crystal. The thick nonlinear crystal works as a frequency filter due to the large group velocity 
mismatch (GVM = l/u 9 (A /2) - l/v g (X Q )) 5]. The working condition for the GRENOUILLE requires that GVM x 
L » r t , where L is the confocal parameter (the length of the crystal can be used for order of magnitude calculations) 
and T t the time duration of the pulse, that can be taken to be the full width at half maximum of the pulse being 
measured. Group velocity dispersion (GVD = l/v g (X — 5X/2) — l/u 9 (Ao + SX/2)) can be avoided by imposing the 
• extra relation GVD x i « t c , where r c is the pulse coherence time (approximately the reciprocal bandwidth, 1/A/) 
These two conditions are combined in Equations (0) and 10, where TBP stands for time bandwidth product. 

'. GVD ] < J < GVM (2) 
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> TBP = — (3) 



Consider a transform limited (TBP ~ 1) input pulse with center wavelength of 800 nm and a FWHM of <5A = 10 
nm. Consider also a 5 mm thick BBO crystal. Equation (J2Jl results in inequality which shows that the crystal is 
appropriate for measuring the input beam [f| . 

GVD = 20/s/cm < 200/s/cm < GVM = 2000 fs/nm (4) 

As a counter example, consider a pulse with TBP ~ 10. Suppose that we have a bandwidth of the order of SX ~ 20 
nm and hence a pulse time length of the order of 100 fs. Equation J2J is not satisfied for the 5 mm BBO crystal. 
Table (1) presents a comparison for the case TBP ~ 10 and center wavelength A c ~ 800 nm. 



'Electronic address: silvaac@physics.umd.edu URL: http://www.wam.umd.edu/~silvaac 
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Spherical lens(f= 100mm) 

A 




Cylindrical lens(f= 100mm) 

Grenouille Side View , Imaging System 

FIG. 1: Grenouille's experimental setup from a side view. The picture is not drawn to scale. Nonlinear crystal: BBO dimensions 
5 x 10 x 10 mm. Fresnel biprism has an apex angle of 168° and almost square base of 1 x 1 in. The camera used is a Pulnix 
TM-72EX with a filter for the fundamental beam adapted to its opening. 
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II. EXPERIMENTAL SETUP 

GRENOUILLE simplest experimental setup is presented in Figure ^ where the vertical axis of a measured FROG 
trace is automatically aligned to be wavelength and the horizontal axis, time delay. Folding mirrors can be included to 
further reduce the final size of the setup (the commercial GRENOUILLE is a small rectangular box of 5cm x 15cm x 
25cm)but that does not change the experimental set in any essential way. 

The alignment has to closely achieve Figure ^ configuration with the second harmonic going in the center of the 
CCD camera and the fundamental beams being blocked by the slit. All the distances shown in Figure^are within 2 
mm for an ideal FROG trace (It seems enough to use a ruler to measure these distances). It is "easy" to get some 
SHG "trace", the problem is to have a good quality trace and also to be sure what has been measured is indeed 
a FROG trace. After the initial alignment, fine adjustments to the imaging system as well as the position of the 
input cylindrical lens have to be done. The distance between the input cylindrical lens and the crystal is found by 
maximizing the intensity of the second harmonic beam after the crystal. The position as well as function of each 
element in GRENOUILLE is now discussed. 



A. Input beam 

The size of the input beam determines the range of wavelength and time delay seen in the FROG trace and 
consequently limits the spectral range as well as time range of the pulse to be measured. 

The spectrometer in traditional FROG apparatus is replaced by the thick non-linear crystal and by the size of the 
input beam. The thick crystal determines the resolution and the input beam the range. 



Side view of non-linear crystal 
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FIG. 2: Neglecting the polarization, this shows that although the crystal was cut for a certain wavelength (say 800nm, the 
center beam with angle a)the beam coming with an angle (3, close to a, will also phase match as long as the pulse has a 
wavelength component that phase matches at (3. If (3 is the most extreme angle in the range of angles of the input beam inside 
the crystal that phase matches, the wavelength range of the setup is the wavelength that frequency doubles between a and 
(3. For this setup the polarization is perpendicular to the surface of the figure for the incoming beam and on the plane of the 
figure for the SHG beams (blue and green beams in figure). 



The range of angles covered by focusing the input beam inside the nonlinear crystal determines the wavelength 
range that can be measured by GRENOUILLE. Figure [21 shows that only a small wavelength range is phase matched 
at each angle (the smallness of wavelength range is determined by the crystal filtering) . This mechanism creates at 
each angle a or f3 a small range of wavelength that is further mapped into the CCD camera. 

Consider the GRENOUILLE constructed with a BBO crystal and designed to work with input beams from 10 mm 
to 20 mm at a center wavelength of 810 nm. Beams with waist of about 10 mm have a range of angles inside a BBO 
crystal of about 1.72°. Phase matching angle for 810 nm is about 28.9°, for 860 nm is about 27.3° From these 

values, the wavelength range with a 10 mm input beam is about 810 ± 25 nm. 

The wavelength range increase of the GRENOUILLE can not be achieved by simply increasing the vertical size of 
the input beam. The confocal parameter has to be considered, since it is the effective length of the crystal. Bigger 
input beam waist also implies smaller confocal parameter. Smaller confocal parameter might also imply violation of 
Expression @ , the fundamental condition for GRENOUILLE to work. 

The delay line in traditional FROG traces is replaced by the horizontal beam waist and the Fresnel biprism. The 
range of the time delay axis is given by the size of the horizontal beam waist crossed inside the crystal through the 
action of the biprism (Fig. |2J). The expression for the delay is just given by the size of the horizontal crossed trace 
inside the crystal (D) and by the angle (9) formed by the crossed beams (Eq. 



At 



D x sin 9 



(5) 



For an input beam size of 10 mm, BBO crystal, and a Fresnel biprism with apex angle 168°, the delay axis is 
At < 2.3 ps. 
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FIG. 3: The effect of the input cylindrical lens and the Fresnel biprism is to create two crossing strips of the input beam inside 
the crystal. The width of these strips is responsible for the time delay in the FROG trace (blue beam). 



The function of the Fresnel biprism [5j is to create from the input beam two beams that cross inside the nonlinear 
crystal (Fig. [3J. The position of the biprism between the focusing lens and the nonlinear crystal can be chosen to 
adjust the range of the delay axis. The maximum delay axis for a given crystal is calculated with Equation [S] where 
the input horizontal beam size plays a fundamental role discussed in the previous section. 

The optimal size of the delay axis does not need to be the maximum delay size achievable with a given 
GRENOUILLE setup. The required condition for the range of the delay axis is to include the majority of the 
FROG trace. The FROG trace has to show all its intensities including its outer edges (horizontal outer edges may be 
clipped for time delay axis that are not wide enough). The smallest delay axis range that shows the full FROG trace 
is the optimum delay axis range. Smaller delay axis give better quality FROG traces since more of the fundamental 
beam is concentrated on the area of the second harmonic conversion increasing the intensity of the trace. 

For an input beam of 10 mm (horizontal waist), the distance between the BBO crystal and the Fresnel biprism 
that maximizes the delay axis is about 100 mm. In practice, one can search for the best position by looking at the 
FROG trace while moving the Fresnel biprism in between the crystal and the focusing lens. 

The Fresnel biprism can be mounted on a fixed optical mount. Rotation on the perpendicular plane of figure ^ or 
tilt towards the nonlinear crystal or the focusing cylindrical lens is not needed, as long as the input focusing cylindrical 
lens can be rotated on the plane perpendicular to the plane of figure ^ The rotation degree of freedom is critical for 
the fine adjustments on GRENOUILLE's alignment. Rotating the input cylindrical lens, helps aligning the reflection 
symmetry about the wavelength axis inherent in the GRENOUILLE FROG traces (SHG FROG trace is even with 
respect to the delay). 



The Nonlinear crystal is chosen to meet the requirements given in Equations @ and |J3}. The combination of the 
input beam and nonlinear crystal is equivalent to a spectrometer where the resolution is determined by the nonlinear 
crystal. 

The second harmonic generated pulse intensity from a nonlinear crystal is given in the frequency space (u>) by 
Equation ©, where L is the length of the crystal, AK the phase matching condition and GVM the group velocity 
mismatch between the fundamental and the SHG beam |9|. The sine 2 (sine = sinx/x) term limits the bandwidth for 
the SHG pulse. The bandwidth becomes narrower with the increase of GVM and also narrower with the length (L) 
increase of the crystal. Hence to increase filtering, maximum GVMxL is a desired quality for the Nonlinear crystal, 
which leads to the fundamental conditions for GRENOUILLE to work (Eqs. (|2I3|) ). 



B. Fresnel Biprism 



C. Nonlinear Crystal 
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(6) 



The bandwidth of the filtering process is calculated using Sellmeier equations 0, Q] and depends on the confocal 
parameter of the focused input beam inside the nonlinear crystal. Most of the SHG takes place along the confocal 
parameter which then defines the effective length of the crystal. The minimum bandwidth in a setup for a given 
nonlinear crystal is achieved with a confocal parameter of the length of the crystal. A BBO crystal 5 mm long with 
the input beam at 800 nm has the minimum bandwidth of about 2.5 nm. The bandwidth for a confocal parameter 
of 2 mm is approximately 2.8 nm even though the crystal is still 5 mm long |5j. 

The most important point in aligning the nonlinear crystal is to set the principal plane on the vertical plane of the 
table (for the set up of Figs. ^3 and [2J> - The polarization of the incoming beam should be parallel to the line focus of the 
incoming cylindrical lens. This will guarantee that the vertical axis is the wavelength axis. It is important to phrase 
this point since one can always have SHG generation, even with the crystal at the wrong orientation. With the crystal 
correctly set, one should see 3 SHG beams at the same height. The FROG trace is the center beam. Translating 
the crystal on the plane perpendicular to the incoming beam might be helpful to avoid defects on the surface of the 
crystal. Tilting and rotating the crystal helps on the alignment since phase matching is angle dependent. 



D. Imaging System 

The function of the imaging system is to take the FROG trace generated by the nonlinear crystal and image it 
on the surface of the CCD array. The combination spherical lens and cylindrical lens (Fig. ^) is assembled back to 
back on a mount with rotational degree of freedom. These combination creates a vertical focus of / = 100 mm and 
a horizontal focus of / = 50 mm. 

The vertical direction of the imaging lens maps the wavelength dependent angle from the FROG trace onto the 
surface of the CCD array as in Figure 0] The position on the camera (x) is a nearly linear relation of the wavelength 
(A). 



CCD 




FIG. 4: Side view of schematic representation of the equivalent lens on vertical axis of the imaging system for the setup of 
figure Q The represented lens is cylindrical and has a focus /=100 mm. The SHG beam coming from the nonlinear crystal 
has wavelength as a function of angle A(il). The imaging lens takes Ai = Const x Af2 = Const x AA. Knowing the center 
wavelength X(Q — 0) from an independent measurement, the wavelength along the x axis can be found. Note that the linear 
relations between x and A are only valid for small angles (Q). 
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The time delay axis is collected by the horizontal direction of the imaging lenses and directed into the CCD camera. 
In this set up (FigU), the size of the FROG trace along the time delay axis on the CCD array is equal to the one 
produced by the crystal. 

The position of the CCD camera is set to have the FROG trace in the center. The FROG algorithm needs as an 
input parameter the center wavelength, which is half of the center wavelength of the pulse being measured only at 
the center of FROG trace. 

The resolution of the Frame grabber and CCD camera has to be higher if the FROG trace is more complex. For 
simple traces (such as Gaussian pulses), having a FROG trace of size 64 x 64 pixels is enough to recover the correct 
phase and intensity of the input pulse [8j]. For this setup (Fig.(@J), the resolution used is 512 x 480 pixels. 

The overall alignment of the GRENOUILLE has to achieve a time delay symmetrical FROG trace located at the 
center of the frame grabber view area. 



The scale for the wavelength axis and time delay axis is found with the aide of an Etalon. The etalon also helps in 
lifting the time direction ambiguity inherent in GRENOUILLE (SHG) FROG An air spaced etalon with spacing 
d produces pulses separated by 2d/c in time, where c is the velocity of light in air. The resultant pulse trace has 
features that can be used to calibrate both time and wavelength axis. The direction of time is also known since the 
second pulse coming out of the etalon is always of lower intensity than the first. Figure [5] shows a FROG trace after 
an air spaced etalon with d — 50 /im. The features used to calibrate GRENOUILLE are the modulations both in 
the horizontal and vertical axis. The modulations on the delay axis are separated by 2d/c. Measuring the number of 
pixels separating the modulations gives the scale for the delay axis (2d/(c x ^pixels)). 

The vertical modulations are related to the horizontal modulations by ArA/ = 1. Using this relation to find the 
wavelength separation of the vertical modulations, one only needs to divide the wavelength separation by the number 
of pixels. 

Typical values for the set up of Figure and for the trace in Figure [5] are, for the horizontal axis (time delay), 8 
fs /pixel and for the vertical axis (wavelength) 0.14 nm/ pixels. 

Errors in the calibration induce errors in the recovered pulse. The size of the error in the recovered electric field is 
at least of the same size of the error in the calibration. An error of 20% in the time delay axis leads to an error of at 
least 20% in the time width of the recovered electric field. The exact influence of the error in the calibration depends 
on the complexity of the pulse to be measured and generally is bigger for more complicated pulses [Tol ] . 



The final goal of FROG is to determine the complex electric field from the FROG trace. This is done with the aid 
of a numerical algorithm that starts with an initial guess of the complex electric field and iterates to the measured 
electric field. The algorithm assumes two constraints. The first is the form of the electric field generated by the 
nonlinearity in use. For SHG, the electric field is given by Equation (J7J). The second constraint is the experimentally 
acquired FROG trace which is given by Eq.(JSJ) or Eq.Q). Extra constrains are reported to create instabilities in the 
algorithm (for instance overflow), including adding an independent measured spectrum as one of the constrains 



The algorithm works as presented in Figure 10- The initial guess can be any electric field, the outcome should 
not depend on the initial guess, but the convergence time is minimized for clever guesses. The commercial code 
(Femtosoft) gives 4 options for initial guesses but does not allow for external input. 

Once the initial guess is chosen the algorithm enters into a loop. The steps in the loop depend on the details of 
the algorithm. There are at least 8 different algorithms that can be used LiUJ The most common can be found in 
the commercial versions of the FROG algorithms (Femtosoft and MarkTech), they are: the Basic FROG (or Vanilla 
FROG, first algorithm created), Generalized Projections, Short-cut Generalized Projections, Projections Over-Step, 
Intensity Constrained Basic FROG, Over-Correction Basic FROG and Multidimensional Minimization Technique 0]. 



E. Calibration 



III. FROG PULSE RETRIEVAL ALGORITHM 



E sig (t,T) =E{t)E{t-r) 



(7) 




(8) 



FIG. 5: FROG trace from a train of pulses separated by multiples of 333.33 fa. Such FROG trace is produced by an air spaced 
etalon in the path of a transform limited 810 nm Spectra-Physics Tsumami Thsapphire oscillator with pulse width of ~ 190 
fs. The horizontal axis is time delay and the vertical axis is wavelength. 

The algorithm that has been reported as most robust and responsible for the success of FROG is Generalized 
Projections Generalized Projection was shown to converge to the real pulse for all FROG traces where the other 
algorithms fail. Generalized Projections works as follows. The first step (step 1), E S i g is constructed with Equation 
Q and then I frog is calculated with a fast Fourier transform (step 2). In step 3 the calculated I frog has its 
magnitude replaced by the magnitude of the measured FROG trace (Eq.(|5J). In step 4, E sig (t,r) is constructed by 
taking the inverse Fourier transform. Until now all the steps are also true for the vanilla FROG algorithm. While the 
Vanilla FROG takes step 5 only by integrating E S i g (t, r) on r to get back E(t), Generalized Projections takes step 5 
by introducing distance minimization (Z) between E sig (t,r) and E(t)E(t — r) to get E(t) back (Ect.QlOjl). Once the 
Step 5 is complete, the cycle is reinitiated. 



E slg (t,r) 



E s jg(t,T) 
\E si g(t,T)\ 



x [Ifrog^^t)} 1 



(9) 




(10) 
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FIG. 6: General diagram for the FROG algorithm adapted from Q . For SHG FROG and GRENOUILLE the signal generation 
is given in Eq.J7J.The loop 1-5 continues until desired error level is reached. 



G =[^ E [lFROG(u,T)-\E sig (u;,T)m /2 (11) 

U.T — 1 

The error in the pulse retrieval (G) is measured with the FROG error (Eq.JTJJ), where N is the grid size of the 
FROG trace and \E S i g (uj, r)| 2 is normalized to a pick value of unity before G is calculated [lflj. The size of the FROG 
error decreases after each loop and typical error values depend on the quality of the trace and on the size of the grid. 
The program converges to the correct pulse for values of the error less than 10~ 4 for noise-free pulses (theoretical 
traces). It is not uncommon to achieve errors less than 0.5% for good quality simple SHG traces 128 x 128 in size. 

The quality of the traces and hence the minimum error (G) level achievable when using the FROG algorithm 
depends on the noise level when recording the trace. For multiplicative noise {Ifrog-Noisc = I frog x (1 + noise)), 
G ~ e x (TBP/N) 1 / 2 , where TBP is time-bandwidth product of the pulse, e is the error in the trace data points 
where the trace is nonzero and N x TV is the array size. For additive noise (I frog— Noise = Ifrog + noise) the error 
isG~e HEI. 

Despite the success of Generalized Projections to find the intensity and phase of a light pulse, this technique, as well 
as any other technique, is not proven to converge for all pulses. The strategy adopted in the commercial programs, is 
to combine several algorithms into the program and to switch between them during a run. The commercial Femtosoft 
algorithm switches between strategy everytime the error does not decrease by 0.5%. This approach hopes to combine 
the strength of all the strategies and eventually cover all the possible FROG traces found in nature. 



A. Practical issues when using Femtosoft FROG pulse recovery program 

If pulse to be recovered is not theoretical, noise filtering is appropriate. It is customary to first remove the lowest 
pixel and then the edge of the data in the window after the trace has been imported into the program. This operation 
should be done only once to avoid clipping the low intensity data. 

After noise filtering the data has to be extracted for the program to initiate its cycle. It is possible to extract only 
a small window around the FROG trace avoiding the background noise by selecting with the mouse the desired area. 
This will result into a smother recovered electric field but the algorithm has a recently discovered bug. This operation 
will rescale the wavelength axis incorrectly. 



9 



Acknowledgments 

I thank Pat O'Shea and Mark Kimmel for the indispensable counsel of how to set up GRENOUILLE. In spe- 
cial, I thank Pat O'Shea for detailed discussion on the basics of FROG and GRENOUILLE. I also thank Bhaskar 
Khubchandani with direct help in setting up GRENOUILLE. Finally, I thank Prof. Rajarshi Roy for advice, guidance 
and logistic support, without which this work would have not been possible. 



[1] D.J. Kane and R. Trebino, Single-Shot Measurement of the Intensity and Phase of an Arbitrary Ultrashort Pulse By Using 
Frequency-Resolved Optical Gating. Opt. Lett., 1993. 18(10): p. 823-825. 

[2] R. Trebino and D.J. Kane, Using Phase Retrieval to Measure the Intensity and Phase of Ultrashort Pulses: Frequency- 
Resolved Optical Gating. J. Opt. Soc. Amer. A, 1993. 10(5): p. 1101-1111. 

[3] Rick Trebino, Frequency-Resolved Optical Gating: The measurement of Ultrashort Laser Pulses, 1st ed. (Kluwer Academic 
Publishers, New York NY, 2002). 

[4] R. Trebino et al., Measuring Ultrashort Laser Pulses in the Time- Frequency Domain Using Frequency- Resolved Optical 
Gating. Review of Scientific Instruments, 1997. 68(9): p. 3277-3295. 

[5] P. O'Shea, M. Kimmel, X. Gu, and R. Trebino, Highly simplified device for ultrashort-pulse measurement. Opt. Lett., 26, 
932 (2001). 

[6] V.G. Dmitriev, G.G. Gurzadyan and D.N. Nikogosyan, Handbook of Nonlinear Optical Crystals, 3rd ed. (Springer, Berlin, 
1999). 

[7] P. O'Shea, M. Kimmel,X. Gu and R. Trebino, Increased bandwidth in ultrashort-pulse measurement using an angle-dithered 
nonlinear-optical crystal. Opt. Expr., 2001. 7(10). 

[8] K.W. DeLong, D.N. Fittinghoff, and R. Trebino, Practical Issues in Ultrashort-Laser-Pulse Measurement Using Frequency- 
Resolved Optical Gating. IEEE J. Quant. Electron., 1996. 32(7): p. 1253-1264. 

[9] J-C. Diels and W. Rudolph, Ultrashort Laser Pulse Phenomena, 1st ed(Academic Press, San Diego CA, 1996). 
[10] K.W. DeLong, et al., Frequency-Resolved Optical Gating With the Use of Second-Harmonic Generation. J. Opt. Soc. Amer. 
B, 1994. 11(11): p. 2206-2215. 

[11] D. J. Kane, Real-time measurement of ultrashort laser pulses using principal component generalized projections, IEEE 
Journal of Selected Topics in Quantum Electronics 4, 278-84 (1998). 

[12] K.W. DeLong, et al., Pulse Retrieval in Frequency-Resolved Optical Gating Based on the Method of Generalized Projec- 
tions. Optics Letters, 1994. 19(24): p. 2152-2154. 

[13] D. N. Fittinghoff, K. W. DeLong, R. Trebino, and C. L. Ladera, Noise sensitivity in frequency-resolved optical-gating 
measurements of ultrashort pulses, Journal of the Optical Society of America B 12, 1955-1967 (1995). 



